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\  ABSTRACT  .  r:  - .  ,  .  : 

A  formulation  for  determining  the  electromagnetic 
field  penetration  through  a  circular  aperture  is  developed 
using  Bab inet’s  principle  and  the  singularity  expansion 
method.  Computational  procedures  for  determining  the 
penetration  field  of  the  aperture  are  discussed  and  a  low 
frequency  check  of  the  procedures  is  proffered. 
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INTRODUCTION 

In  the  study  of  various  problems  of  electromagnetic  compatibility 

it  is  necessary  to  determine  the  penetration  field  of  an  aperture.  At 

low  frequencies  a  quasi  static  solution  is  available,1  and  at  high 

2 

frequencies  the  Kirchhoff's  solution  is  available.  Also  a  rigorous 

3 

solution  developed  by  Flammer  is  available.  He  applies  Babinet's 

principle  and  determines  the  scattering  from  the  disk  equivalent  of  an 

aperture  in  a  conducting  plate  by  considering  the  disk  as  a  limiting 

form  of  an  oblate  spheroid  and  constructs  vector  wave  function  solutions 

to  the  Helmholtz  wave  equation  in  oblate  spheroidal  coordinates.  However 

no  numerical  results  are  presented;  but  it  is  expected  that  the  greatest 

amount  of  energy  penetrates  the  aperture  near  the  first  resonances  of  the 

aperture,  the  intermediate  frequency  region. 

Recently  the  singularity  expansion  technique  has  been  found  to  be 

,  4 

useful  in  the  study  of  electromagnetic  pulse  interaction.  In  principle 
one  should  be  able  to  take  Flammer's  solution  and  develop  singularity 

4 

expansions  for  the  induced  current  and  charge, as  Baum  did  in  treating 

the  sphere.  '  However,  identifying  and  calculating  the  natural  frequencies 

and  modes  may  be  difficult.  An  alternative  treatment  used  here  is  to 

develop  an  integral  equation  for  the  induced  current  density  on  a  Babinet 

equivalent  disk,  solve  the  integral  equation  numerically  using  the  method 

of  moments,  and  apply  the  singularity  expansion  method  numerically  similar 
5 

to  Tesche. 


2 


ANALYSIS 


7Z  I  '-'tVj  ‘1 


--’•■The  solution  for -the  electromagnetic  field  penetrating  an  aperture 

in  a  conducting  plane  may  be  determined  by  solving  the  complimentary  disk 

3 

problem  and  applying  Babinet's  principle.  The  electromagnetic  form  of 
Babinet's  principle  states  that  if  (Ej^H^)  is  the  scattered  field  when 
(Ej,Hj)  is  the  field  incident  in  the  positive  z  direction  on  a  perfectly 
conducting  disk  lying  in  the  plane  z=0,  then 


=  -  C  Sf  ,  H®  -  i-  E®  ,  z  -  0 

E®  -  x,  H®  ,  H®  =  -  i  E®  z  -  0 

are  the  diffracted  fields  when  the  wave 


*2 


4  wf 


(la) 


(lb) 


is  incident  in  the  positive  z  direction  on  the  complimentary  perfectly 

conducting  screen  with  an  aperture.  Here  C  - Vy/e  is  the  intrinsic  wave 

Impedance  of  the  medium  surrounding  the  aperture  or  disk.  In  the  case  of 

the  aperture  the  total  field  in  the  half  space  z  -  0  is  formed  from  the 

superposition  of  the  incident  wave,  the  reflected  wave  in  the  absence  of 

an  aperture,  and  the  diffracted  field  (E®,H|).  In  the  half  space  z  -  0, 

is  the  total  field  penetrating  the  aperture. 

To  determine  the  scattered  or  diffracted  field  from  the  aperture  it  is 

necessary  to  obtain  the  induced  surface  current  density  on  the  disk.  The 

Integral  equation  for  the  surface  current  density,  3g,  on  a  disk  centered  at 
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the  origin  and  in  the  z»0  plane  is 
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t;  •  •  -  ^2  ypjf 


+  grad '  /  div  J8 ( r) G (r | r  * ) dS 


(2) 


A  A 


A  A 


where  tj  “  r'- and  ■  ♦'  with 


G(r  |r') 


|r-r'| 


The  incident  field  may  be  expanded  in  transverse  magnetic  (TM)  and 
transverse  electric  (TE)  cylindrical  modes.  For  the  TM  case  there  is  no 

i  ' 

magnetic  field  component  along  the  z-axis  and  for  the  TE  case  there  is  no 
electric  field  component  along  the  z-axis.  An  arbitrary  electromagnetic 
field  may  be  expressed  as  the  linear  combination  of  a  TM  part  and  a  TE 
part  (see  Appendix  1).  These  two  cases  are  treated  separately. 

TM  Mode  Excitation 

( 

The  form  of  the  TM  mode  expansion  for  the  incident  plane  wave  field  is 


"V  :  l 
,  "l  S  -» 


E 

m-o 


E®  (r)  cos  m<j>  r 

L  l,n 


-  i  '  -i 

+  E°  (r).sin  m$  *  e"Jkz  cos  el  (3) 

Z,m  *  J 


where  6^  is  the  angle  of  incidence,  the  angle  between  the  direction  of  pro¬ 
pagation  and  the  direction  of  the  positive  z-axis.  The  corresponding  form 
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of  the  induced  current  distribution  on  the  disk  is 


J,(r)  -  £ 

m*0 


[k®  ' (r)  cos  m$  r  +  K°  (r)  sin  m$  $ 

L  l,o  z,m  J 


(4) 


-tv, 


< 


Using  the  representations  (3)  and  (4)  in  (2)  yields  coupled  integral 
equations  for  K|  ffl(r)  and  m(r)  to  be  solved  for  each  value  of  m.  They  are 


E!..<r’>  ■  i 


”  l/[Vl<r|r)  +  G"fl<t|r’>]  K>."(r)rdr 

a 

+  7  Sr  f 


-  W'l'’)]  K2,m 


(r)rdr 


and 


2m  d  r  ,  n 

+  3  7T  f  G_(r  r*)K° 

dr*  /  m  1  2,m 


(r)dr 


(5) 


E°  (r’)  = 

2,m  4tt 


j*y*[Gm-l(r!r,)  ~  Gnri-l(rlr’)]  Kt,m(r)rdr 


?•  / 


Gm(r|r’)(ri-+  1)K®  (r)dr 
dr  xtm 


a 

+  /'[G»- 1  (r  |r ■  1  (r  I  r ■ )]  K° ^ (r)  tdr 


2m 

k 


2  a 


f-,  f  Gm(r|r')K° 

2r*  I  °  2,m 


(r)dr 


(6) 
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where 


/  -jk  V(r-r')2+2rr'(l-cos<j 
Jm(rjr')  =  /  - -  —  - - - 

v- y  ^  (r-r')2+2rr’  (l-cos$)* 


■  :  7 ,  <. 

(7) 


Apparently  Gm  as  shown  above  possesses  a  logarithmic  singularity  at  r-r\ 
Since  (7)  must  be  evaluated  numerically  a  more  convenient  form  is  needed. 
After  some  mathematical  manipulation  (7)  becomes 


3m(r|r’)  =  2^  j[« 


e'JkK  cos(2m0)-l  /R  (  d0 


where 


|(r-r')2  +  4rr’  sin2  0 


and  K  is  the  complete  elliptic  integral  of  the  first  kind.  Note  that  the 
integral  in  (8)  does  not  possess  a  singularity  at  r=*r’  and  may  be  evaluated 


numerically./,- 


3 r/sns. 


A  different  singularity  appears  to  occur  when  r'»0  and  r=0  unless  a 
restriction  is  placed  on  the  current  expansion.  The  evaluation  of  (7)  at 
r'-O  yields 


-  e-jkr  - 

G»(r l0)  "  *  6mo  (9) 

where  6mo  13  the  usual  Kronecker  delta.  In  order  for  the  respective  integrals 
in  (5)  and  (6)  to  remain  integrable  at  r-0  the  following  is  required 


K?,o<r) 


r  +  fl 


(const. )r 
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The  physical  equivalent  to  (10)  Is  the  requirement  of  a  finite  charge 
distribution  at  r=0. 

Another  restriction  to  be  satisfied  by  the  current  distribution  is 
the  radial  component  of  the  current  must  vanish  at  the  edge  of  the  disk. 

That  is, 

Ke  (a)  =0  (11) 

1 ,  m 

for  all  m.  Further  note  that  one  may  set 

K°  Q(r)  -  0  (12) 

Because  of  the  singularities  in  Gm  the  derivatives  of  the  integrals 
must  be  evaluated  carefully.  If  the  derivatives  are  taken  inside  the  integrals 
the  integration  is  performed  by  using  Cauchy  principal  values.  To  avoid 
this  additional  complication  the  derivatives  of  the  integrals  will  be 
evaluated  numerically  by  using  finite  differences. 

TE  Mode  Excitation 

6 

The  form  of  the  TE  mode  expansion  for  the  incident  plane  wave  field  is 

::  ••  .  -  ;  ~.j.»  r  '  \  ^ '  •'  ■  i'.h  A 

itf(r)  *  I  E°  (r)  sin  m$  r 

.  3  m=0  L  1,0  •  ---v .-3  -  .  V  ; 

-  E®  (r)  cos  + 1  e"Jkz  co8  9i  (13) 

z»m  J 

6 

And  the  corresponding  form  of  the  induced  current  distribution  on  the  disk  is 

3„(r)  =  Z  |k°  (r)  sin  m$r  -  K5  (r)  cos  m$  $  (14) 

m=0L  l*m  Z,m  J 
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Using  the  representations  (13)  and  (14)  in  (2)  yields  coupled  integral 

equations  for  K°  (r)  and  K®  (r)  to  be  solved  for  each  value  of  m.  The 
l  ,m  L  ,m 

form  of  the  integral  equations  is  exactly  the  same  as  (5)  and  (6)  when  the 
following  transformations  are  made: 


Ke  (r) 

— (r) 

1  ,m 

l,m 

K °  (r) 

- >K®  (r) 

2,m 

2,m 

E®  (r) 

/-s 

M 

0« 

T 

l,m 

l,m 

E°  (r). 

— >E®  (r) 

2,m 

2,m 

(15) 


In  order  to  have  a  finite  charge  distribution  at  r=0 


K®  (r) - (const. )r 

2‘°  r  -*-  0 

Another  restriction  to  be  satisfied  by  the  current  expansion  i3 


(16) 


1,« 


(a) 


(17) 


For  convenience  one  may  set  -  • 

K^0(r)  =0  (18) 

'-j  ■  "  ' 

for  all  r. 
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NUMERICAL  SOLUTION 


Because  of  the  inherent  complexity  of  the  integral  equations  for 

the  induced  surface  current  density  an  analytic  solution  appears  to  be 

virtually  impossible  to  obtain.  Therefore  a  numerical  solution  technique, 

7 

namely  the  method  of  moments,  is  employed.  Moment  expansions  that  satisfy 


the  restrictions  on  the  induced  surface  current  densities  are 


Ki,„« 


vrrP 


iKrjrp+l.rp) 


If  (19)  and  (20)  are  substituted  into  (5)  and  (6),  and  the  resulting 
equations  are  satisfied  at  r=rp,  p=l,2,...N,  a  system  of  linear  equations 
are  obtained  for  the  ctp's  and  the  Bp's.  This  procedure  is  sometimes  referred 
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to  as  point  matching  or  colocation.  The  resulting  system  of  equations 


( 


is 


2N+1 

Z  n(m,I,J)F(m,J)  «*  r(m,I) 
J*1 


where  I-1,N  and  J*1,N 


n(m,I,J)  =-  —  I  Fl(m-l,I,J)  +  Fl(nrt-l,I,  J) 

Ar  L 

+  F2(m-l,I,J)  +  F2 (m+l,I,J) 


!  ^F5(m,I+l, 


J)  -  F5(m,I,J) 


k(Ar)4 

+F6 (m, 1+1, J)  -  F6(m,I,J)J 

+rJ_1[F3(m,I+l,J)  -  F3(m,I,J)] 

[FACm.I+l.J)  -  F4(m,I,J)]J 


J+l 


for  1*1, N  and  J=1,N+1 

H(m, I,  J+N)  =»  -  I  FI (m-1 ,1,  J)  -  Fl(m+l,I,J) 

L 

+  F2 (m-l,I,J)  -  F2(nrH, 

,J)  -  F3(m,I,J) 


F3(m,I+l,: 


+  — — j 
k(Ar) 

+  F4(m,I+l,J)  -  F4(m 
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(21) 


for  1=1, N+l  and  J=1,N 


Il(m,  I+N,  J) 


Fl(m-l,I,J)  -  Fl(mH,I,J) 

Ar  L 

+  F2(m-1 ,1, J)  -  F2 (m+l,I, J) 

—  |  2|F5(m,I,J)  +  F6 (m, I ,  J) 
kAr  (  L 


+  rj_j  F3(m,I,J) 


for  1=1, N+l,  J=1,N+1 


-  r 


J+l 


F4(m,I,J) 


] 


i 

I 


and 


II(m,  I+N,  J+N)  = 


4 


Ar 


2m 

kAr 


Fl(m-l,I,J)  +  FI  (nrt-1 ,1,J) 

+  F2(m-l,I,J)  +  F2(m+l,I,J)J 

j^F3(m,I,J)  +  F4(m,I,J)j 


F(m, J) 


(m) 

aJ 

e(m) 

J-N 


J=l  ,N 
J=N+1,2N+1 


(22) 


r(m,I) 


E2,m^rI-N^ 


1=1, N 

I=N+1,2N+1 


(23) 


The  integral  functions  FI,  F2,  F3,  F4,  F5  and  F6  are  defined  in  Appendix  II. 

It  should  be  observed  that  the  foregoing  applies  for  both  TM  and  TE  excitation. 
The  appropriate  modal  expansion  for  the  incident  fields  are  given  in  Appendix  I. 
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At  this  point  the  induced  surface  current  density  may  be  obtained  by 
solving  the  foregoing  system  of  linear  equations  using  the  digital  computer. 
Obviously  a  system  of  equations  must  be  solved  for  each  modal  current. 
Andreasen6  who  treats  the  body  of  revolution  suggests  that  the  maximum 
number  of  modes  needed  is 

m|  ‘v  ka  sin  8^  +  6  (24) 

max 

for  ka  sin  8^  £  3  and  m|max  1  as  0^  0. 
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SINGULARITY  EXPANSION  METHOD 


According  to  the  singularity  expansion  method  the  natural  frequencies 
may  be  obtained  by  searching  for  the  zeros  of  the  determinant  of  the  fore¬ 
going  system  matrix  JI(m,k).  From  this  point  forward  in  the  analysis  the 
Laplace  transform  frequency  variable  s  •  jck  is  used  for  the  frequency. 
Thus  the  natural  frequencies  are  obtained  from 


det 


II(m,s) 


-  0 


(25) 


for  frequencies  sa  independent  of  the  index  m.  Note  that  the  natural 
frequencies  for  the  TE  modes  are  the  same  as  for  the  TM  modes. 

The  solution  for  the  induced  current  distribution  is 


-1. 


F(m,s)  =  II  (m,s)r(m,s)  (26) 

for  each  mode.  Applying  the  singularity  expansion  method  (26)  becomes9* 

M„(b)  t 


F0»,s)  -  E  C‘(B)r(m,s) 

aS”S  01 

rv 


(27) 


where  M  (m)  is  defined  as  the  natural  mode  vector  and  is  the  solution  to 


the  equation 


n(m,sa)Ma(m)  =  0 


(28) 


and  Ca(m)  is  referred  to  as  the  coupling  vector,  and  satisfies  the  equation 


II  (m,s  )C  (m)  =  0 

a  a 


(29) 


The  M  (m)  and  C  (m)  are  normalized  according  to1* 


c>> 


fen(m,s) 


M„(m) 


(30) 


s=*sr 


Here  the  class  II  form  for  the  coupling  coefficient  is  used. 


Actually  F(m,s)  represents  the  mth  cylindrical  mode  contribution  to 
the  current  distribution..  Summing  the  cylindrical  mode  contributions, 
see  (4)  and  (14),  yields  the  final  form  for  the  current  distribution.  The 
foregoing  (25)-(30)  apply  for  both  TE  and  TM  excitation.  For  TE  excitation 
use  the  TE  mode  field  components  in  T  as  defined  in  (13)  and  (23)  and  for 
TM  excitation  use  the  TM  mode  field  components  in  T  as  defined  in  (3)  and 
(23). 

To  obtain  the  time  domain  response  of  the  Induced  current  distribution 

5 

on  the  disk  the  appropriate  Laplace  transform  of  (27)  is  evaluated.  It  is 


F(m, t)  -  Z  U(t)M  (m)C*(m)r(m,aa)l 
a 


where  U(t)  is  a  diagonal  square  matrix  of  unit  Heaviside  functions  which 


serves  to  enforce  the  requirements  of  causality. 
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SCATTERED  FIELD 


Once  the  induced  current  distribution  of  the  disk  is  known  then  the 

8  . 

scattered  field  may  be  readily  determined.  The  electric  field  is 

-  j  kR~ 

E® (r* )  =  -  j_JL_  ^..7)7  +  k2?g]  | -  ds  (32) 


and  the  magnetic  field 


)ds 


(33) 


where- 


-jkR  - 


(34) 


and 


.  -jkR  ~  i  , 

(Js*7)7  | - -  -k2(Js-R)R  +  |-(jk+|) 


[3<JS-WR-Is]|i 


,-jkR 


(35) 


with  R=r'-r  and  R=R/R.  Note  that  Jg  is  given  by  (4)  or  (14)  depending  upon 
the  mode  of  excitation.  The  far  field  approximations  to  (32)  and  (33)  may 
be  readily  determined  by  neglecting  the  R”  and  R  J  terms  appearing  in  (34) 
and  (35). 

In  order  to  verify  the  numerical  procedures  that  are  employed,  low 
frequency  excitation  may  be  considered  and  the  equivalent  dipole  moments  of 
the  scattered  field  obtained.  These  are 


P  -  —  f 


JgdS 


(36) 
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(37) 


Mo  =  7  J  r  x  Jsds 


s 

From  using  a  quasi-static  approximation  it  is  found  that1 

PQ  «  a3  jf(x*Ednc)x  +  (y*Einc)y J  (38) 

Mq  -  -  A  a3(z‘Hinc)z  (39) 

substituting  the  current  representations  (4)  and  (14)  into  (36)  and  (37) 
yields 

Pox  "  j£/[Kl.l(r>  -  K2,l(r)]  tdr  (40> 

o 

<4l) 

Moz  =  -  */0  K2,0(r)rdr  (42) 


At  low  frequency  (40)-(42)  may  be  compared  with  (38)  and  (39)  to  check 
the  numerical  solutions  for  the  components  of  Jg. 

The  field  penetrating  an  aperture  in  an  infinite  plate  may  be  determined 
by  applying  the  electromagnetic  form  of  Babinet's  principle.  As  discussed 
previously  the  field  transformations  (la)  and  (lb)  are  used  with  the 
scattered  field  from  the  Babinet  equivalent  disk  to  determine  the  field 
diffracted  by  an  aperture.  The  scattered  field  from  the  disk  is  given  by 
(32)  and  (33). 
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CONCLUSION 


The  problem  of  the  electromagnetic  field  penetration  through  an 
aperture  In  a  perfectly  conducting  plate  la  formulated  by  using  Bablnet's 
principle.  A  numerical  solution  for  the  Induced  current  distribution  on 
the  Bablnet  equivalent  disk  Is  developed.  This  numerical  solution  Is 
formulated  both  as  a  direct  moment  method  solution  and  a  singularity 
expansion  solution.  Finally  computational  procedures  for  determining 
the  penetration  field  of  the  aperture  are  discussed  and  a  low  frequency 
check  on  the  computational  procedures  is  proffered. 
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APPENDIX  I 


For  convenience  consider  a  plane  wave  propagating  in  the  direction 
0a  H  -  and  to  be  polarized  in  the  direction  forming  an  angle  0p  with 
the  unit  vector  -0.  Here  0  is  the  usual  polar  angle  of  the  spherical 
coordinate  system.  Therefore  the  incident  electric  field  may  be  resolved 
into  two  components 


where 


E*  =  -  E|'  0  +  E^  y 

_ ii  „  .  -ik(z  cos  0.-x  sin  0.) 

E  =  E_  cos  0_  e  J  '  i  i' 

1  o  P 

Ef  =  E  sin  0  e-Jk(z  cos  6i~x  sin  6i> 


(Al) 


(A2) 


(A3) 


Then  by  the  principle  of  superposition  each  component  of  J1  is  considered 
to  be  a  plane  wave.  The  magnetic  field  associated  E**  will  he  e  no  z 

x  1 

component  -  a  TM  wave.  And  the  z  component  of  E^  will  be  zero  -  a  TE  wave. 


TM  case 

If  the  electric  field  E^  is  expanded  in  a  Fourier  series  and  set 
equal  zero,  then  (Al)  yields  - 

00 

-  E 

m**0 

+  E°  ffl(r)  sin  o$  e  COS  01 

where 


D 


E‘,»« 


cos  m$  r 
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-  cos  0  cos  0 .  el 
p  i 


E1  m(r) 
l,m 


j'(kr  sin  0.)  E 
m  i'  o 


E2  m(r) 

2  fin 


cos  0  cos  6,  e  1 
p  i  nr 


nH-1  Jm(kr  sin  V 
m  — : - - — - 


kr  sin  0, 


with  (r,<j>,z)  as  the  usual  cylindrical  coordinates  and  e  as  Neumann' 

m 

number. 


TE  case 

If  the  electric  field  is  expanded  in  a  Fourier  series  and  e| 
equal  zero,  then  (Al)  yields 


E  I  E°  (r)  sin  m<l>  r 
m=0  L  1*“ 


+  E®  (r)  cos  m<|>(()-j  e  ^kz  cos  6i 
2,m  J 


where 


O  nu.1  J  (kr  sin  0, ) 

E1  =  “  sin  6n  1,n  - Z - 5 - A -  E 

l.m  p  kr  sin  0^  o 


E2,m<r> 


n>4*1  » 

sin  0  eJ  x  J  (kr  sin  0. )  E 
p  m  m  I  c 
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s 


I 

set 


ft  -  *ca  i  -Trffr'Tjtflih  I 


APPENDIX  II 


The  various  Integral  functions  used  in  H(m,I,J)  are  defined  here 
Fl(m.I.J) 

For  J  1+1,  I>1 


where 


Fl(m 


Ar 

,I,J)  -  J  ( 


(u+rJ_1)Gm(rI lu+rj_i)udu 


Ar 


+2 


r  illioLK 

J  ri+rj-i+u 


rI~rJ- L~u 
rI+rJ-l+u 


du 


V»lv>  ■  2  J  j[« 


cos(2m9) 


-']/R| 


de 


o 

R  »  (u-v)  +  4uv  sin^G 
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[  2  4  6  8 1 
ao  +  alc  +  a2C  +  a3?  +  a4?  J 

r  2  4  6  8 "I.  2 

-|b0  +  b1?Z  +  b2?  +b3C  +b4C  J*nC 


rx  =■  (I-l)Ar 
rj  •  (J-l)Ar 
Ar  »  a/N 


a3 


1.386 

294 

4 

bo  "  °*5 

0.096 

663 

443 

bx  =  0.124 

985 

94 

0.035 

900 

924 

b2  =  0.068 

802 

486 

0.037 

425 

637 

b3  -  0.033 

283 

553 

0.014 

511 

962 

b4  -  0.004 

417 

870 
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otherwise 


i,I, J)  =  J 


<r„T-u)G«(riirj+ru)ud“ 


/1"S,4r  (rJtl-»)u  [~rI~rJ+l+u 

J  rI+rJ+ru  rI+r"l'" 

r>  u  “* 


For  I»1 


F2(m,I, J)  = 


+  (Ar)Z  qx 


ijkAr(l-1kAr)-l1  e“jkrJ+l 
k* 


otherwise 


For  J  j*  I  and  J  ^  1+1 


rir 

-  V 

O 


I lu+rJ-i)du 


»r— s— « u 

/  ri+rj-i+u  L 

»  r> 


ri-rj-i-u 

ri+rj-i+u 


For  J-I+l 


F3(m, 


GmCrI  i^j-] 
J  o 


rI+rJ-l+u  I  J-l 


rI~rJ-l~u 
r_+r_  .+u 


For  J-I 


F3(m 


rAr 

,I,J)  =  VrjIufr^Jdu 


+  2/ 


,( 1-6 )Ar 


rI+rj_l+u 


ri-fj-ru 

ri+rj-i+u 


Ar 

+  ^  ,i 


F4(m,I,J) 

For  J  I  and  J  4  1+1 


rAr 

F4(m,I,J)  -  I  Gm(ri  !rj+i“u)du 
J  o 

+  2  /‘Ar  1 

/  ri+rj+ru 

»  n 


ti'rj+i+" 


ri+rj+i-» 


du 


For  J-I 


F4(m,I,J)-l  Gm(rI|rJ+1-u)du 


r- 

rAr 

J  6Ar 


ri+rj«-" 


ri~rj+i+u 

rt+rj+ru 


du 


^  Ar 
+ 
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For  J-I+l 


Ar 


F4(m,I,J)  =  f  G  (r  |r 

/  m  I '  J+l 


-u)du 


(l-6)Ar 


+  2  r  — i — K 

/  ri+rj+ru  1 

ft  1— 


rI"rJ+l+u 


rI+rJ+l“u 


du 


+  — — ■  q 
ri  1 


FSQn.I.J) 

For  J  4  I,  I>1 


F5(m,I,J)  = 


/Gm(ri  lrJ_1+u)udu 
o 


+2J 

j  o 


Ar 


ri+rj-i+u 


K 


rI+rJ-l+u 


udu 


For  J=I,  I>1 


r 

F5(m,I,J)  =1  Gm(rIjrj  ^Judu 
J  o 

r 


I-5)Ar 


+  2 

. "  a  ... 


ri+rj-i+u 


ri~rj~i~u 

rI+rJ-l+u 


udu 


+  ^  ,, 


C 
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